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ABSTRACT 



Context. The environmental conditions for asympotic giant branch (AGB) stars to reach the carbon-rich (C-rich) phase 
are important to understand the evolutionary process of AGB stars. The difference between the spatial distributions 
of C-rich and oxygen-rich (O-rich) AGB stars is essential for the study of the Galactic structure and the chemical 
evolution of the interstellar medium (ISM). 

Aims. We quantitatively investigate the spatial distributions of C-rich and O-rich AGB stars in our Galaxy. We discuss 
the difference between them and its origin. 

Methods. We classify a large number of AGB stars newly detected by the AKARI mid-infrared all-sky survey. In the 
color-color diagrams, we define their occupation zones based on the locations of known objects. We then obtain the 
spatial distributions of C-rich and O-rich AGB stars, assuming that they have the same luminosity for a given mass-loss 
rate. 

Results. We find that O-rich AGB stars are concentrated toward the Galactic center and that the density decreases 
with Galactocentric distance, whereas C-rich AGB stars show a relatively uniform distribution within about 8kpc of 
Sun. 

Conclusions. Our result confirms the trends reported in previous studies and extends them to a Galactic scale. We 
discuss the relations between our result, the Galactic metallicity gradient, and the chemical evolution of the ISM in our 
Galaxy. 

Key words. Galaxy: structure - Stars: AGB - ISM: chemical evolution 
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1. Introduction 
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Asymptotic giant branch (AGB) stars are low- to 
intermediate-mass (1-8 M©) stars at their final evolution - 
ary stages that undergo mass- loss fe.g. lHabing et al.lll985lL 
Their ages range from 0.1 to several Gyr, depending on 
their initial masses. The spatial distribution of AGB stars 
roughly traces past star-formation activities. 

The AGB stars are divided into three classes based on 
their carbon-to-oxygen elemental abundance ratios: carbon- 
rich AGB stars (C-rich; C/0>1), oxygen-rich AGB stars 
(O-rich; C/0<1), and S stars (C/ Q~l). Many AGB stars 
in the Milky Way are O-rich (e.g. lHabing et al.l [l985t) and 
show silicate-dominated circumstellar features. In contrast, 
some fraction of intermediate-mass AGB stars (1.1-4 M ) 
become C-rich, displaying the features of both carbon- 
bearing molecules (e.g. C2, CN, HCN, C2H2, and CO) and 
carbonaceous dust (e.g. amorphous carbon and graphite). 
This chemical change is triggered when carbon atoms, 
which are synthesized in the stellar core, are dredged up 
to the surface by the deep convective envelope. 



The difference in the spatial distributions between O- 
rich and C-rich AGB stars, i.e., the variation in the ra- 
tio of C-rich AGB stars to O-rich AGB stars (C /M ra- 
tio) f rom place to plac e, has long been studied (|Blancol 
119651: IWesterlundl Il965l ). It is known that O-rich AGB 
stars concentrate toward the Galactic center (G.C.), 
while C-rich AGB stars are distributed more uniformly 



around Sun (e.g. Thronson et al.l 1987 ]_ Claussen et al.l 



Il987t Uura fc Kleinmannl Il990t iNoguchi et al.l I2004T ). For 



example 
tra of 689 
rections (/ 



Le Bertre et al.l ( 2003lT investigated the spec- 
acquired in the four radial di- 
130°, and 170°) with the Near- 



AGB stars 
=-10°, 50°, 

Infrared Spectrometer (NIRS; iNoda et aL 119961) on-board 
the I nfraRed Telescope in Space (IRTS: iMurakami et al.l 
119961 ). and found that O-rich AGB stars outnumber C-rich 
ones in the inner Galaxy (<8kpc), whereas the situation is 
the reverse in the outer Galaxy. However, most past studies 
surveyed limited regions of the solar neighborhood because 
of their sensitivity constraints. A comprehensive study of O- 
rich and C-rich AGB stars covering the entire Milky Way 
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and a quantitative discussion of their dependence on envi- 
ronment (e.g. metallicity) has not been performed yet. 

The galactic AGB stars are believed to be a significant 
contributor to the interstellar medium (ISM). Their distri- 
bution and amount of mass loss are interesting in view of 
the chemical evolution of the ISM in our Galaxy and the 
life cycle of heavy elements in the universe. Recently, su- 
pernova remnant (SNR) is focused as the candidate to fill 
the missing budget i n the supply o f the interstellar dust 
(e.g. iRho et al.ll200l iLee et al-lbool llshihara et aTll2010ft . 
In particular, many discussions have focused on the supply 
of dust from the remnants of core-collapse (Type II) su- 
ernovae since the d iscovery of dust in the early universe 
Bertoldi et al.ll2003h . However, the supply of carbonaceous 
grains should be discussed independently of the supply of 
the total amount of dust grains dominated by silicates to 
clarify the whole picture of the life cycle of ISM and explore 
clues to the origin of life. The gas and dust budgets and 
their impact on the ISM evolution were quantitatively in- 
vestigated bv lMatsuura et all (|2009T ) using sensitive Spitzer 
images of the Large Magellanic Cloud (LMC), and the total 
mass loss amount from C-rich AGB stars revised. 

From 2006 to 2007, AKARI (jMurakami et al.ll2007[) sur- 
veyed the entire sky with six infrared photometric bands 
between 9-160 (tm. The mid-infrared (MIR) part of the 
all-sky survey was perf ormed with two bro ad bands cen- 
tered at 9 and 18 /im (jlshihara et al.l l2006f ) using one of 
the on-boar d inst ruments, the InfraRed Camera (IRC; 
lOnaka et al.l I2007D . The detection limits (5cr) for point 
sources per scan are 50 and 90mJy for the 9 and 18 ^m 
bands, respectively, with spatial resolutions of about 5", 
thus surpassing the IRAS survey in the 12 and 25 (im bands 
by an order of magnitude in both sensitivity and spatial res- 
olution. More than 90% of the entire sky was observed at 
both bands during the lifetime of the liquid helium cryo- 
gen. In total, 870,973 sources (844,649 in the 9 /jm band and 
194,551 in the 18 /jm band) were extracted with the uniform 
thresh old and compiled as t he first p o int-so urce catalogue 
(PSC; llshihara et aTll2010l ). Ilta et al.l (|2010[) reported that 
the MIR catalogue contains at least 90,000 AGB stars. 

In this paper, we investigate the spatial distributions 
of C-rich and O-rich AGB stars across the Galaxy. We de- 
scribe the data analysis applied to select C-rich and O- 
rich AGB samples from the AKARI/MIR all-sky survey 
PSC (hereafter AKARI/MIR PSC) and the reliability of 
the samples in Sect. 2. The spatial distributions of and 
any differences between C-rich and O-rich AGB stars in 
our Galaxy are presented in Sect. 3. Section 4 discusses 
the origin of the difference in the distributions of C-rich 
and O-rich AGB stars, compares our results with related 
and extragalactic results, and describes the effects on the 
chemical evolution of the ISM in our Galaxy. A summary 
is given in Sect. 5. 



2. Data analysis 

Source classifications based on the loci of sources in color- 
color diagrams were often applied to previous sets of in- 
frarcd survey 'data for IRAS (e.g. IWalker fcCohenl 119881: 



irared survey data 
IWalker et all [19891 



(e-g- 
(e-g- 



van der Veen fc Habingl 11988^ MSX 



an et all 120011: iLumsden et al.l 1200 



Buchanan et al 



and Spitzer 
method basi- 



l2009h . This type of 
cally utilizes the characteristics of the spectral energy dis 
tributions (SEDs) of the detected sources. It is an effee 



tive and unique approach in classifying a large amount of 
unidentified objects with only photometric data, though 
contamination is inevitable. We apply this method to the 
AKARI/MIR PSC, which contains many unidentified ob- 
jects, and classify all the detected objects, selecting C-rich 
and O-rich AGB stars among them. 

We initially cross-identify all 9 fim and 18 /im sources 
in th e AKARI/MIR PSC f31 data set (jlshihara et al.l 
I2010D with the 2MAS S Point Source Catalogue version 7 
(jSkrutskie et al.|[2006l ) based simply on given celestial coor- 
dinates. The search radius of 3" was optimize d for the cross- 
ident ification between these two catalogues (jlshihara et al.l 
I2010D . In total, 614,204 AKARI/MIR sources are cross- 
identified with the 2MASS catalogue. We also search 
for the identification of the AKARI/MIR sources in the 
SIMBAD database. The search radius is 10". In total, 
331,764 AKARI/MIR sources have plausible identifica- 
tions. In other words, more than 60% of the AKARI/MIR 
sources are new compared to the SIMBAD database. 

All the cross-identified sources arc plotted on color-color 
diagrams. Figures QJi, [lb, and [TJ; show the color-color di- 
agrams of [J]-[K] versus (vs.) [K]-[9], [K]-[9] vs. [9]~[18], 
and [J]-[K] vs. [9]-[18], respectively. Here, [J], [K], [9], and 
[18] represent Vega magnitudes in the 2MASS J, 2MASS 
K s , AKARI 9/im, and AKARI 18 /im bands, respectively. 
The zero point magnitudes for th e 9/im and 18 urn bands 
are 56.3 and 12.0 Jy, respectively (jOnaka et al.ll2007j) . 

We then define eight groups (A-H) and assign already 
known objects to the groups following the identification of 
the SIMBAD database as summarized in Tabled] Figures^ 
[3[ and[4]show the distributions of the known objects on the 
same color-color diagrams as Figs. Hb, and [Tb, respec- 
tively. In all three diagrams, a large fraction of objects are 
concentrated near the origin, which represents the colors of 
the photospheres of normal stars. The object distribution 
extends to the upper-right direction for stars with warm 
circumstellar dust. In the diagrams including the [J] — [K] 
axis (i.e. Figs. []Ji and Q}:), [K]-excess objects are mixed 
with objects reddened by interstellar extinction. In the di- 
agrams including the [K] — [9] axis (i.e. Figs. [TJi and [TJd), 
very red ([K] — [9]>5) objects, presumably dominated by 
hot dust emission or the emission of polycyclic aromatic hy- 
drocarbons (PAHs), can be clearly distinguished from the 
majority of the objects, although they are relatively small 
in number. 

We define the occupation zone of each class exclusively 
on these color-color diagrams based on the locations of 
the known members. The defined zones are overlaid in 
Figs. El El andSJ Table [1] summarizes the numbers of known 
members used in the definition of the occupation zones, 
and the total number of the sources in the corresponding 
zones. The occupation zones defined in this work are con- 
sistent with the criteria used in the source selecti ons in the 
related works using AK ARI/MIR survey data (jlta et al.l 
I20T0I: lTaku~a et al.ll2010h . 

The C-rich and O-rich AGB stars can be clearly distin- 
guished in all the three color-color diagrams (Figs. [2^, Ek, 
and Hk). This is because the silicate features of O-rich 
AGB stars around 10 and 20 /im are well-covered by the 
AKARI 9 fim and 18 /im bands, while the carbonaceous 
(e.g. SiC) features of C-rich AGB stars do not contribute 
to these AKARI bands (Fig. [5|. These characteristics are 
effective in helping us to differentiate C-rich from O-rich 
AGB stars, and their availability is one of the unique advan- 
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tages of the AKARI/MIR survey data compared with th e 



other infrared surveys o f IRAS (iNeueebauer et al 



MSX (|Price et all l200l . and WISE dWright et al 



1984) 



2010) 



The distinction of C-rich from O-rich objects is clearest in 
the [J]-[K] vs. [9] -[18] color-color diagram (Fig.SK). Thus, 
we first use the classification in the [J] — [K] vs. [9] — [18] 
color-color diagram (hereafter called candidate samples), 
and then treat the sources residing in the occupation zones 
in all the three color-color diagrams as more purified sub- 
samples (hereafter called purified samples). As a result, the 
candidate samples contain 18,596 C-rich and 83,012 O-rich 
AGB stars, while the purified samples contain 5,537 C-rich 
and 11,416 O-rich AGB stars (Table 0). 
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Fig. 5. Relative response curves of the AKARI 9{im 
band and 18/im band. The ISO/SWS spectra of typi- 
cal Gal actic C-rich and O -rich stars are overlaid as ref- 
erences (jSloan et al.ll2003j) . An example of the C-rich star, 
IRC +40540, displays carbonaceous features around 11 and 
30 fim, that are not covered by the AKARI bands. A typi- 
cal O-rich star, oCet, shows silicate features around 10 and 
18 /im, which arc efficiently covered by the AKARI bands. 



As described above, contamination is inevitable for this 
type of method. We investigate the reliability and the com- 
pleteness of the classified samples as follows and summarize 
them in Table The samples located in each occupation 
zone are composed of four types of objects: (a) correctly 
classified known objects; (b) known objects that are clas- 
sified as other categories, but can be regarded as the same 
group; (c) unknown objects; and (d) known objects that 
should belong to another group. The objects labeled as 
variable stars and infrared objects are included in (b) for 
C-rich and O-rich stars. Mira-type stars are also included 
in (b) for O-rich stars. For each occupation zone, we de- 
fine (n a + Tib)/(n a + rib + n c + n<i) as a lower limit, and 
(n a + rib + n c )/(n a + rib + n c + rid) as an upper limit to 
the reliability of the classified samples, where n a , rib, rz c , 
and rid are the numbers of objects in (a), (b), (c), and (d), 
respectively. For each object class, the completeness is es- 
timated as the ratio of the number of correctly classified 
known objects to the total number of known objects of the 
same class in the AKARI MIR PSC. The number of known 
objects for each class depends on the previous incomplete 
surveys. However, that is only the method to assess our 
samples, that is useful to judge the validity of the statis- 



tical results based on our samples. For the candidate sam- 
ples, the reliability is 44-95% for C-rich AGB stars, and 
53-95% for O-rich AGB stars. The completeness is 82% for 
C-rich AGB stars and 80% for O-rich AGB stars. For the 
purified samples, the reliability is 68-96% for C-rich AGB 
stars, and 71-96% for O-rich AGB stars. The completeness 
is 64% for C-rich AGB stars and 29% for O-rich AGB stars. 
The completeness for O-rich AGB stars is lower than for C- 
rich AGB stars. As shown in Figures [21 El and SI since the 
O-rich and C-rich samples arc bluer, they are more inter- 
mixed, whereas the redder ones in the same samples can 
be clearly separated. Therefore, with the present classifica- 
tion, it is inevitable that the completeness for the redder 
members is higher than that for the bluer members. Thus, 
the completeness of the C-rich samples is higher than that 
of the O-rich samples, which contain a larger numbers of 
lower-mass non-dusty stars. 

3. Results 

3.1. Galactic distributions 

Figure [5] shows the spatial distributions of the objects se- 
lected in the [9] — [18] vs. [J] — [K] color-color diagrams in 
density maps in the Galactic coordinates. The surface den- 
sities in the maps in Fig. [5] are calculated for the area of 
lxl deg 2 , and given in units of deg -2 . The spatial dis- 
tributions of S-stars, C-rich AGB stars, and O-rich AGB 
stars trace the Galactic disk. The bulge, the LMC, and the 
Small Magellanic Cloud (SMC) arc also recognized in these 
maps. Dwarfs and giants (group C in Fig. [5]) show relatively 
uniform distributions in contrast to the other populations 
except for two dim concentrations around I = —90° and 
+90°. Since normal stars are faint in the MIR bands, their 
detection volume is limited to the local space (< lkpc) 
around Sun. Th e distributions are sim ilar to that of the 
Hipparcos stars (jPerrvman et al.lll997t ). corresponding to 
the structure of the local arm. In the panel for YSOs (group 
A in Fig. [6]), nearby star- forming regions such as Orion, 
Taurus, and pOph arc recognized. These maps demonstrate 
the plausibility of the classification by reproducing well- 
known major characteristics of the spatial distribution of 
each Galactic object. 

3.2. Distributions of C-rich and O-rich AGB stars 

Figures [7^ and [T]b show the celestial distributions of the 
purified C-rich and O-rich AGB samples, respectively, as 
density maps in Galactic coordinates. By assuming that 
the star counts can be described using Poisson statistics, 
the 2-er significance levels of the surface densities of the 
C-rich and O-rich AGB samples are 0.3 and 1.5 deg -2 , re- 
spectively. The O-rich stars show a concentration toward 
the G.C. (the areal density varies from 3 deg -2 in the G.C. 
to <0.5deg -2 in the outer Galaxy (|Z| > 90°)), whereas C- 
rich stars show a uniform distribution along the Galactic 
plane (the areal density is kept to within 0.2-0.4 deg -2 
across the Galactic plane). Thus, the decrease in density 
along the Galactic plane seen in Fig. [TJd is significant, while 
the small-scale structures in Fig. [?K are not. These results 
can be accounted for either by an intrinsic difference in the 
distribution on a Galactic scale, or by a situation where O- 
rich stars are distributed on a Galactic scale in the thin disk 
with interstellar extinction and C-rich stars arc located in 
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relatively nearby regions. To investigate any difference in 
greated detail, three-dimensional information would be in- 
dispensable. 

We plot C-rich and O-rich AGB stars in the projec- 
tion onto the Galactic plane of both disk stars < 10°; 
Figs. |Sh. and c) and stars at high Galactic latitudes (\b\ > 
10°; Figs. [SJd and d). The distances of the individual sam- 
ples were estimated from their AKARI 9 /im fluxes (-Fg^m) 
and [K] — [9] colors. A formula to estimate the distances was 
empirically derived using the sub-samples for which the dis- 
tance and the mass-loss rates were known. These parame- 
ters were estimated based on MI R observations indepen - 
dently by IZhang et all (j2010t ) and iLe Bertre etall (j2003h . 
We derived the relation among the distance, the mass-loss 
rates, and Fgum for the sub-samples. We also derived the 
relation between the mass-loss rates and [K] — [9] colors (see 
Appendix for details of the derivation) . By combining both 
relations, we obtained a formula to derive the distance from 
Fg^u! and [K] — [9]. The maximum distance that we can 
probe in our sample is limited by the detection limit for the 
AKARI 18 /im flux. The maximum distance to the faintest 
members in our sample is estimated to be 8±2 kpc, based on 
the above flux-distance relation with typical [9] — [18] colors 
(also see Appendix). 

The edge of the outer part of our Galaxy (the left side of 
Fig. [SJ indicates the end of the corresponding stellar popu- 
lations, whereas the edge toward the direction of the oppo- 
site side over the G.C. (the right side of Fig. [8]) corresponds 
to the detection limit of this survey for the corresponding 
objects. The elliptical concentration is seen around the G.C. 
in the panels for the disk component (|fe| < 10°; Figs. [8^, 
and [5b). The inclination of about five degr ees to the Sun- 
G.C. axis may indicat e the central bar ([Genhardl l2002t 
ICole fc Weinberell2002D . However, we note that the elon- 
gated distribution is also caused by the uncertainties in the 
derived distances to the individual stars. The concentra- 
tions of sources around I ~ 280-300° in the panels for high 
latitude objects (|6| > 10°; Figs. [5b and[Sbl) indicate the 
LMC and the SMC. They are clearly separated from our 
Galaxy, though there remains an elongated distribution to- 
wards the LMC presumably caused by the errors in the 
distance estimates. 

The most remarkable finding is that O-rich AGB stars 
are concentrated toward the G.C. while C-rich AGB stars 
have a relatively uniform distribution across the Galactic 
disk. 

4. Discussion 

To have more quantitative discussion about the difference 
in the spatial distributions of C-rich and O-rich AGB stars, 
we plot the arcal number density of our samples as a func- 
tion of Galactocentric distance (R g ) in Fig. [9l At various 
galactocentric distances, we count the numbers of disk AGB 
stars included in the annular regions of width 0.3 kpc within 
8 kpc of the Sun considering the survey depth. The smaller 
errors indicated by darkly colored stripes in Fig. [S] are esti- 
mated by the Poisson statistics in the star counts, while the 
larger ones indicated by more lightly colored stripes include 
an uncertainty of ±2 kpc in estimating the distance. We also 
investigate the completeness of both samples as a function 
of R g in the lower panel of Fig. [9] The completeness of the 
C-rich sample increases with R g , while that of the O-rich 
sample is lower in the outer galaxy (R g > 9 kpc) . For O-rich 



AGB stars, the completeness is higher for the redder mem- 
bers (see Sect. 2). The lower completeness in outer Galaxy 
means that the fraction of dusty members (i.e. massive O- 
rich AGB stars) becomes smaller in the outer Galaxy. In 
contrast, the completeness of the C-rich AGB stars shows 
little depends on their color. The number and variety of 
the catalogued C-rich AGB stars are relatively small in the 
outer Galaxy. Thus, the completeness of C-rich AGB stars 
increases toward the outer Galaxy. The ratio of the com- 
pleteness of the C-rich sample to that of the O-rich sample 
is found to be almost constant at R g < 9 kpc. Hence, our 
result is not affected much by the variations in the com- 
pleteness. 

Table [3J summarizes the related investigations with 
the numbers and the penetrating depths of the samples. 
Previous studies have consistently reported relatively flat 
distributions of C-rich AGB samples compared with O-rich 
AGB samples, which are instead concentrated toward the 
G.C. Our result in Fig.[9]quantitatively confirms the trends 
reported in the previous works, which improves the spatial 
completeness by extending the survey volume to 8 kpc from 
the Sun and avoiding source confusions toward the G.C. 




J) I I I I I I I MlllJ 

2 A fl & 1D 12 U 



Fig. 9. Upper panel: areal densities of the number of the pu- 
rified C-rich (red) and the O-rich (blue) AGB samples in the 
Galactic plane (\b\ < 10°) as a function of Galactocentric 
distance (R g ), not corrected for the completeness (see be- 
low). For each R g , we count the numbers of the disk AGB 
stars included in the annular regions of width 0.3 kpc within 
8 kpc of Sun (the penetrating depth of our sample) and di- 
vide the star counts by the corresponding areas. The nar- 
rower stripes indicate the errors estimated from the Poisson 
statistics for the star counts, while the wider stripes indi- 
cate the ones derived by assuming uncertainties of ±2 kpc 
in estimating the distance. Lower panel: completeness of 
both samples as a function of the Galactocentric distance 

Why do the distributions of C-rich and O-rich AGB 
stars differ ? One possible reason is that the metallicity 
of the ISM at their birth place affects the stellar chemical 
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composition of the photospheres (e.g. lNoguchi et aT1l2004[) . 
The metallicity gradient along R g was reported by many 
authors (e.g. by using Cepheids and open clusters). The 
reported gradient of d[Fc/H]/dR g at R q = 5-15 kpc varies 
from -0.04 to -0.13 dcx kpc" 1 (Pedicelli et al.ll2009D . and 
is higher toward the inner Galaxy. The [O/H] gradients are 
also r eported by several authors (e.g. Smartt &; Rollestor] 
Il997t IShaver et al.l Il983l iRudolph et al.ll2006t) . to be be- 
tween about —0.04 and —0.07 dex kpc -1 at radii 6-18 kpc. 
The log(C/M) value at radii R g =2 to 14 kpc in our result 
has the gradient of ~0.1dex kpc -1 , which has an inverse 
gradient of slope similar in absolute value to that of the 
metallicity. 

The anti-correlation between the metall icity and the 
C /M r atio is also reported for nearby galaxies. lBrewer et al.l 
(|l995h reported on M31, a spiral galaxy similar to our 
Galaxy, observing five regions located at 4-32 kpc from the 
G.C. The metallicity decreases along the R a , whereas the 
C/M ratio increases. iPritchet et al.l ( 19871 ) confirmed that 
the number of C-rich AGB stars per unit luminosity shows 
a correlation with the average [Fe/H] among the galaxies, 
IC1613, SMC, NGC6822, NGC55, M33, LMC, NGC300, 
M31, and our Galaxy, thus the lower metallicity results 
in the higher C/M ratio. The relatio n is approximated by 

rf(fFe/Hl)/rf(loK (C /M)) 0.5 (e.g. iBattinelli fc Demers 1 

l2005tlCioni 1120091) . Our result in Fig.|Ushows that the slope 

is d([Fe/H])/d(log(C/M)) = -0.4 1.3 at R=2-14kpc, 

which connects smoothly with the value for nearby galaxies. 
Thus, the trend reported in nearby galaxies is applicable to 
our Galaxy, while the ranges of metallicity are different: 
-1.3<[Fe/H]<0 for nearby galaxies, and -0.5<[Fe/H]<l 
for our Galaxy. We note that the metallicity dependence 
of the C/M ratio in nearb y galaxies is based on optical 
(IBattinelli fc Demers 112005]) or near- infrared (NIR; [C lorn 
[2009F observations, whereas our result is based on MIR ob- 
servations. Figure|H]also suggests that C/M~ 1 at [Fe/H]=0 
(i.e. the solar abundance), although our O-rich samples are 
biased toward dusty members. 

How docs the difference in spatial distributions between 
C-rich and O-rich AGB stars affect the ISM in our Galaxy 
? Our result indicates that as dust suppliers the popula- 
tion of C-rich stars are minor contributors in the inner 
Galaxy but major contributors around the Sun and in the 
outer Galaxy. It is reasonable to assume that the mass loss 
from AGB stars has a direct effect on the chemical com- 
position of their local ISM. The spatial constraint on the 
dominant replenishment of carbonaceous grains to the ISM 
may have implications for habitable zones in our Galaxy 
(jGuillermo et al.ll200lD . though a significant theoretical gap 
remains between the supply of carbonaceous grains to the 
interstellar space and the origin of organic material. The 
spatial variation in the amount of interstellar carbonaceous 
grains relative to that of silicate grains has not been re- 
ported on the Galactic scale. Furthermore, it is well-known 
that the PAH (as a representative of carbonaceous grain) 
emission displays a good spatial corre l ation with the cold 
dust emission (e.g . lOnaka et al.l Il996l iDraine et al.l l2007t 
iBendo et al.ll2008T ); the latter is mostly emitted by silicate 
grains. We plan to perform detailed and quantitative com- 
parisons of the spatial distributions of carbonaceous and 
silicate grains with those of C-rich and O-rich AGB stars 
in our Galaxy in future work. 



5. Summary 

On the bassis of AKARI observations, we have investi- 
gated the spatial distributions of C-rich and O-rich AGB 
stars in our Galaxy. The C-rich and O-rich AGB stars in 
the AKARI 9 ^tm and 18 /jm All-Sky Survey catalogue are 
well-distinguished by the empirical approach using color- 
color diagrams. The reliability and the completeness of the 
sample selection are discussed. The C-rich AGB stars are 
uniformly distributed within 8 kpc from Sun, while O-rich 
AGB stars are concentrated toward the G.C. Our result 
confirms those of previous investigations, extends the sur- 
vey volume to the Galactic scale (~8 kpc), and enhances the 
completeness toward the inner Galaxy by spatially resolv- 
ing individual sources. We also propose that the metallicity 
gradient with galactocentric distance can be one of the pos- 
sible origins of the difference in the spatial distributions of 
C-rich and O-rich AGB stars. 
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(a) [J]-[K] vs. [K]-[9] (b) [K]-[9] vs. [9]-[18] (c) [J]-[K] vs. [9]-[18] 




Q 2 4 D 10 -10 12 0^0 -1 012345 



t«H8] BH131 [9]-[18] 

Fig. 1. Color-color diagrams in density maps for the AKARI/MIR sources that have 2MASS measurements: (a) [J]-[K] 
versus [K]-[9], (b) [K]-[9] versus [9]-[18], and (c) [J]-[K] versus [91— [18 ] color-color diagrams. The a rrow on each panel 
shows the interstellar extinction vector for Ay = 20 mag. by using the IWeineartner fc Draind (|200ll ) Milky Way model 
of Ry = 3.1. The color scales of the maps are logarithmic in number per unit color space. (The color versions of this and 
subsequent figures are available in the online journal.) 




2 4 6 8 0246 80 246 80 2468 

[K]-[9] [K]-[9] [K]-[9] [K]-[9] 

Fig. 2. Distributions of the known members (identified in the SIMBAD database) shown in contours on the [J]-[K] versus 
[K]-[9] color-color diagram : (a) C-rich AGBs (red; group E), AGBs (black) and OH/IR stars (blue) (O-rich AGBs; group 
F); (b) T Tauri (brown) and YSOs (green) (group A); (c) Be stars (including Wolf-Rayet stars; brown) and emission-line 
stars (green) (group B); (d) normal stars including O-M giants and dwarfs (green; group C); (e) S-stars (green; group 
D), and post- AGBs and planetary nebulae (brown; group G); and (f) normal and starburst galaxies (green) and Seyferts 
(brown) (group H). The occupation zones for the groups A-H are summarized in (g). The contours show 20%, 50%, and 
80% of the peak value. The SIMBAD definitions of the members of the groups A-H are summarized in Table [TJ 
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Fig. 3. Same as Fig. [3J but the data are instead presented on the [K] — [9] vs. [9] — [18] color-color diagram. The occupation 
zone for emission-line stars (B) is not defined on this diagram because they cannot be clearly separated from the other 
objects (panel c). 



i 



(a) 




A 




«* 


l<f° / 

L^L\\( o / 






)<?/o 
\ \ 


w/ 








v o \ 








w% \ 




i 


i 


1 \ 1 


1 




2 3 
[9]-[18] 



012340123401234 
[9]-[18] [9]-[18] [9]-[18] 



Fig. 4. Same as Fig. [2j but shown on the [J] — [K] vs. [9] — [18] color-color diagram. 
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Table 1. Classification of the AKARI/MIR catalogue sources. 



Group 


Representative population 






Number of members 








(Definition in SIMBAD" ) 


[J]-[K] vs 
known 6 


[K]-[9] 


[K]-[9] vs. 
known 6 


[9]-[18] 


[J]-[K] vs. 
known 6 


[9]-[18] 






total c 


total c 


total c 


A 


YSOs (YSO, TT*, HII) 


1,302 


81,802 


776 


17,369 


776 


5,588 


B 


Emission line stars (Be*,Em*,WR*) 


1,984 


1,644 


744 


d 


744 


1,493 


C 


Normal stars (O-M dwarfs & giants; *) 


167,829 


288,597 


35,117 


36,906 


35,117 


27,836 


D 


S-stars (S*) 


852 


205,255 


461 


17,454 


461 


20,973 


E 


C-rich AGBs (C*) 


5,197 


202,252 


3,037 


13,944 


3,037 


18,596 


F 


O-rich AGBs (OH*, AB*) 


1,193 


30,711 


745 


71,190 


745 


83,012 


G 


Post- AGBs & PNe (pA*, PN) 


836 


6,858 


670 


1,762 


670 


2,076 


H 


Galaxies (G, EmG, H2G, rG, LSB, 
LIN, Sy2, Syl, AGN, SyG) 


1,323 


2,735 


744 


3,977 


744 


3,028 



Notes. W Abbreviated expressions for object types used in the SIMBAD database. We assign ambiguous objects to each class. 
For example, the objects labeled as C*? (carbon star candidate) are included in group E. ' 6 ' The number of known objects in 
the SIMBAD database. ^ The number of objects in the occupation zone defined in this work. ' d > Emission-line stars cannot be 
separated clearly on the [K]-[9] vs. [9]-[18] diagram. 



Table 2. Reliability and completeness of the newly classified C-rich and O-rich samples. 



Selection C-rich AGBs O-rich AGBs 

Number Reliability" Completeness 6 Number Reliability" Completeness 6 



(%) (%) (%) (%) 

[J] - [K\ vs. [K\ - [9] 202,253 WM 78 307712 48=92 30" 

[K] - [9] vs. [9] - [18] 13,945 55-89 85 71,191 52-94 71 

[J] - [K] vs. [9] - [18] c 18,597 44-95 82 83,013 53-95 80 
All the diagrams'* 5,537 68-96 64 11,416 71-96 29 



Notes. ( "> The upper limit is calculated for the extreme case that the classification of all the unknown objects located in each 
occupation zone is true, and the lower limit is calculated for the extreme case that the classification of all the unknown objects in 
each occupation zone is false. The ratio of the number of correctly classified known objects to the total number of the known 
objects in each of C-rich and O-rich AGBs. Candidate samples. Purified samples residing in the occupation zones in all the 
three color-color diagrams (see Sect. 2). 



(Group A) YSOs 



(Group C) Normal stars 



(Group D) S-stars 




Fig. 6. Spatial distributions of the AKARI/MIR sources in density maps in the Galactic coordinates of the Aitoff 
projection. These sources are classified on the basis of the [9] — [18] vs. [J] — [K] color-color diagram alone (candidate 
samples). We show groups A, C, D, E, and F because they have a sufficiently large number of members. The color scales 
are linear from 0% to 50% of the peak value and given in units of deg -2 . The coordinate grids and the locations of the 
Galactic bulge, the LMC, the SMC, and the nearby star-forming regions are shown in the lower right panel. 
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(a) C-rich AGBs (purified samples) (b) O-rich AGBs (purified samples) 




Fig. 7. Spatial distributions of the purified samples of (a) C-rich AGB stars and (b) O-rich AGB stars in density maps in 
the Galactic coordinates. The purified sample corresponds to those belonging to all the occupation zones of the [J] — [K] 
vs. [K]-[9] (Fig. [K]-[9] vs. [9] — [18] (Fig.Ot), and [J]-[K] vs. [9] -[18] diagrams (Fig. UK). The color scales are linear 
from 0% to 100% of the peak value and given in units of dcg~ 2 . 



Table 3. Related investigations of the spatial distribution of C-rich AGBs. 



References 


Number 


Distance 


Sample selection 




of sample 


(kpc) 




Clausscn ct al. (1987) 


215 


1.5 


TMSS photometries and optical spectra 


Thronson et al. (1987) 


619 


5 


IRAS colors 


Jura & Kleinmann (1990) 


126 


2.5 


Optical, TMSS, and IRAS colors 


NoEuchi et al. (2004) 


2,250 a 


3.6 


IRAS colors 


This work (2010) 


18,597 i ' 


8 


AKARI/MIR colors 




5,537 c 


8 


AKARI/MIR colors 



Notes. Contamination is estimated to be 40%. ' ' Candidate sample. Purified sample. 
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(a) C-rich ACBs al low Galactic latitudes (|/;| < 10°) (b) G-rich AGRs at high Galactic latitudes (\b\ > 10°) 




-5 5 10 15 20 -10 10 20 30 40 

X (kpc) X (kpc) 



(c) O-rich AGBs at low Galactic latitudes (|6| < 10°) (d) O-rich AGBs at high Galactic latitudes (\b\ > 10°) 




-5 5 10 15 20 -10 10 20 30 40 

X (kpc) X (kpc) 

Fig. 8. Spatial distributions of the C-rich and O-rich AGB stars projected onto the Galactic plane viewed from the north 
pole: (a) C-rich AGBs at |6| < 10°, (b) C-rich AGBs at |6| > 10°, (c) O-rich AGBs at \b\ < 10°, and (d) O-rich AGBs 
at \b\ > 10°. The red dots indicate the purified samples. The contours in panels (a) and (c) show the areal densities for 
the purified samples (i.e. densities of red points). The levels of the contours are 10, 20, and 30 kpc~ 2 for panel (a), and 
10, 50, and 100 kpc -2 for panel (c). The distance scale for each object is estimated from the AKARI 9 /jm fluxes and 
[K] — [9] colors (see Sect. 3). The Sun is located at the origin. The G.C. is to the right of the Sun and shown by the cross 
at (X, Y) = (8.5, 0). The dotted squares in (b) and (d) indicate the size of (a) and (c). Note that the scale of the distance 
used in panels (b) and (d) is a factor of two larger than that in panels (a) and (c). 
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Appendix A: Distance estimate 

We first derive an empirical relation between the distance 
(D; kpc), the observed 9 /xm flux (Fg,,, ™ ; Jy) , and the [K] — [9] 
color for well-stud ied AGB samples dLe Bertre et al.ll2003t 
IZhang et a,l.ll2010h . We then apply this relation to our sam- 
ples and estimate D from the observed values of Fg^ m and 
[K] — [9] . We use AKARI 9 /an fluxes because they are more 
acurate than the 18 (jm f luxes . 

The iLe Bertre et al.l (|2003l ) sample contains 126 C-rich 
star s and 563 O-rich st ars observed by the IRTS/NIRS. 
The IZhang et al.l ([2010D sample contains 184 C-rich stars 
and 110 O-rich stars selected from the literature. Most 
members of both samples have AKARI 9 /im fluxes. The 
mass-loss rates for both samples are estimated observa- 
tionally by using methods independent of MIR observa- 
tions, while the distances are estimated based on NIR 
and MIR fluxes. The mixture of these samples (here- 
after BZ sample) comprehensively covers both dust-poor 
members (10 -7 5 — 10~ 55 M Q yr~ 1 ) and dust-rich members 
(1Q-5-5 _ 10- 4 M Q yr- 1 ). 

Figure IA.1I shows the distances vs. the AKARI 9 (im 
fluxes for the BZ sample. The observed flux is inversely 
proportional to the square of the distance, 



F, 



9fim 



oc L 



9/zrn 



x D~ 



(A.l) 



where L 9lim is the intrinsic luminosity of a star at 9 fim. We 
assume that -Lg^m is a function of the amount of circum- 
stellar dust, thus a function of the mass-loss rate (dM/dt; 
Moyr" 1 ) as 



L 9fim cx (dM/dty, 



(A.2) 



where 7 is a free parameter. 

By combining Eqs. (A.l) and (A.2), the distance is de- 
scribed as a function of F 9fim and dM/dt as 

log D = d+C 2 \og(dM/dt) - i log(F 9Mm ), (A.3) 

where C\ and C2 are obtained as 6.1±0.06, and 0.35±0.01, 
respectively, by fitting the linear relation among log(D), 
log (dM/dt), and log(F 9 ^ m ) for the BZ sample. 

To obtain dM/ dt from our observations, we assume that 
it is a function of excess emission indicated by the [K] — [9] 
color. Figure IA.2I shows the mass-loss rate plotted against 
[K] — [9] color for the BZ sample. We assume that the dotted 
curve in Fig. lA.2[ is an empirical relation between the mass 
loss rate and the color expressed as 

\og(dM/dt) = log(3.8 x (([A'] - [9]) - 0.4)) - 8.0. (A.4) 

From equations (A.3) and (A.4) with observed F 9lim and 
[K] — [9], we estimate the distances of our sample. The re- 
sults are overlaid with dots in Fig. IA.ll The faintest F 9lim of 
our sample is limited by the detection limit of the AKARI 
18 fim band. Typical [9] — [18] colors for dust-poor AGBs 
are 0.5±0.5mag from Fig. [3^. By applying them to the flux 
detection limit at 18 /im (90 mJv: llshihara et al.ll2~010l ). the 
minimum F 9lim for our sample is 270t l^o m Jy, w hich corre- 
sponds to the distance of 8±2 kpc from Fig. IA.ll Figure |A~T1 
also indicates that the AKARI/MIR survey can probe dust- 
rich AGB stars beyond 40 kpc, covering the objects in the 
LMC and the SMC. 



o 




Zhang+(10) O-rich 
r &Bertre+(03) O-rich 
Zhang+(10) C-rich 
LeBertre+(03) C-rich 



6 8 
[K]-[9] 



10 12 14 



Fig. A.2. Mass-loss rate of AG B stars as a functi on of 
the [K] — [9] col o r, me asured by IZhang et al.l (|2010f ) and 
ILe Bertre etaTI (|2003l ). The symbols are the same as in 
Fig.lATl 
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(a) C-rich AGB stars 
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(b) O-rich AGB stars 
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Fig. A.l. (a) Heliocentric distance of C-rich AGB stars as a function of the AKA RI 9 [im flux . The open squares 
and the open ci r cles s how previously measured distances for well-studied samples from lZhang et all (|2010l squares) and 
iLe Bertre et al.l (|2003l . circles), respectively. The black dots show distances estimated for our purified samples. The dashed 
curves indicate the luminosity distance for the objects of [K] — [9] = 14 mag and [K] — [9]=1 mag, which correspond to the 
dust-rich case (dM/dt ~ 10~ 4 M^yr™ 1 ) and the dust-poor case (dM/dt ~ 10~ 8 M Q yr _1 ), respectively (see Fig. \A.2\i . 
The vertical lines indicate the faint flux end of our sample (270^qq mJy), which corresponds to the distance of 8±2kpc 
for the dust-poor members, (b) The same plot as panel (a), but for O-rich AGB stars. 
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